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We present well-controlled results on the metal to insulator 
transition (MIT) within the paramagnetic solution of the dy- 
namical cluster approximation (DCA) in the two-dimensional 
Hubbard model at half-filling. In the strong coupling regime, 
a local picture describes the properties of the model; there is 
a large charge gap A « U. In the weak-coupling regime, we 
find a symbiosis of short-range antiferromagnetic correlations 
and moment formation cause a gap to open at finite tempera- 
ture as in one dimension. Hence, this excludes the mechanism 
of the MIT proposed by Slater long ago. 



Introduction In this letter we report the study of the 
metal-insulator transition (MIT) and its relation to an- 
tiferromagnetism (AFM) in the two-dimensional (2D) 
Hubbard model at half-filling. Results from numeri- 
cal simulations QJ^| have convincingly shown that the 
ground state of the model is an AF insulator, with the 
Neel temperature constrained to be T/v = by the 
Mermin- Wagner theorem. However, the nature of the 
MIT is less clear, as there are two conflicting opinions 
concerning the MIT and its relation to AFM. The first 
opinion is that the strong and weak coupling MI 

transitions are very different. When the local Coulomb 
repulsion parameter U is larger than the non-interacting 
band width W, the ground state is an insulator with a 
large charge gap A w U. The MIT occurs well before the 
onset of magnetism at the temperature T g w U, and the 
spin and charge degrees of freedom are decoupled. The 
superexchange interaction couples the spins with the ex- 
change constant \J\ ~ At 2 /U , and the spins govern the 
low-energy physics. This type of MIT, which is purely 
due to local correlations, is called a Mott transition. In 
weak coupling, a spin density wave (SDW) instability 
develops at T — because of the nesting of the Fermi 
surface. The MIT is the direct consequence of the Bril- 
louin zone folding generated by magnetic ordering. This 
type of MIT will be referred to as a Slater transition. It 
is believed that this regime can be well described by the 
usual many-body weak-coupling approaches. The sec- 
ond opinion, is due to Anderson J(| who has argued that 
the strong-coupling behavior presented above occurs for 
both strong and weak coupling so that a Mott gap is 
present for all U > as in one dimension. As the tem- 
perature falls, local moments develop first because of the 
MIT and then they order so that AFM is the consequence 
of the MIT, not the converse. Because of this strong in- 
teraction, there is no adiabatic continuity between the 



non-interacting and interacting eigenstates. The con- 
ventional renormalizable many-body perturbation theory 
cannot describe this physics. 

The nature of the MIT has been addressed in the Hub- 
bard model in infinite dimensions 0-0] by applying the 
dynamical mean-field approximation (DMFA) The 
unfrustrated model is an example of the first scenario 
presented above with a Slater AF insulator in weak cou- 
pling, and a Mott MIT at strong coupling. The one- 
dimensional model is an example of the second scenario, 
since there is a gap in the charge excitation spectrum for 
all non-zero U p], but no long range magnetic order. 

In 2D, the Mermin- Wagner theorem precludes any 
AFM transition at any finite temperature. Hence, if there 
is a MIT at finite T, it cannot be attributed to the Slater 
mechanism. Finite size simulations (FSS) employing the 
quantum Monte Carlo (QMC) technique have been ap- 
plied to study this model quite extensively during the last 
decade 0-|]l|]. But the difficulty with FSS is that at low 
temperatures, the correlation length is greater than the 
lattice size. Thus the effects of correlations are overes- 
timated for smaller clusters because they are artificially 
closer to criticality than a system in the thermodynamic 
limit. This tendency is only reduced by increasing the 
cluster size, which moves the system in the direction of 
the thermodynamic limit. The spurious AFM gap that 
opens at finite T in the FSS renders the disentanglement 
of the MI and the AF transitions practically impossi- 
ble. The situation is radically different in the dynamical 
cluster approximation (DCA), which is an extension of 
the DMFA. In the DCA, the system is already in the 
thermodynamic limit; however, the DCA restricts corre- 
lations to within the cluster length. As the cluster size 
increases, longer range correlations are progressively in- 
cluded. The FSS and DCA results are complimentary 
in that the DCA (FSS) always underestimates (overesti- 
mates) the effects of non-local correlations |]f3| , |l4j| . Thus, 
the DCA will always systematically underestimate the 
gap formation inherent in the Anderson picture of the 
MIT. Therefore, if a gap is found in the DCA for a finite 
cluster, it will persist in the thermodynamic limit. 

Using the DCA, we show that a MIT occurs in the 
2D Hubbard model at half-filling and finite temperature 
at both weak and strong coupling regimes. This result 
is to be contrasted with the weak-coupling approaches 
which predict that a gap will only exist at T = as a 
consequence of AFM. We argue that at weak couplings 
the symbiosis of short-range AFM correlations and local 
moment formation is the key mechanism for the opening 
of a gap. The local moments forms at a relatively high 
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T leading to a short-range AFM order, that in turn en- 
hances the moments which enhance the AF order, etc. 
This leads to the destruction of the Fermi liquid state 
found in the DMFA studies O]. We discuss the valid- 
ity of our conclusions drawn from finite clusters at the 
thermodynamic limit. 

Method A detailed discussion of the DCA formalism 
was given in previous publications |]l6| , |l7| , ^4[ . The com- 
pact part of the free energy is coarse-grained in reciprocal 
space, projecting the problem onto a finite-sized cluster 
of N c points embedded in a self-consistently determined 
host. When N c = 1, the DCA recovers the DMFA, and 
when N c — > oo, it becomes exact. The cluster problem is 
solved with a generalized version of the Hirsch-Fye QMC 
algorithm |Is|| , and the spectra are analytically contin- 
ued with the maximum entropy method fll9[ ] . We use the 
DCA to study the single band Hubbard model in a two- 
dimensional square lattice at half-filling. The model is 
characterized by W and U . We set W = 2 and we vary 
U from U = W/4 to U = 2W. 
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FIG. 1. The double occupancy D for U = 1 and U = 3 
(inset) for various N c . 

Results Jarrcll , Georges and Krauth Q and Rozen- 
berg et al. || originally obtained a numerically exact so- 
lution from the DMFA in the limit of infinite dimensions 
using QMC. The analysis of these mean- field equations 
revealed the existence of a transition between a para- 
magnetic state and an AFM state. The paramagnetic 
solution was shown to be a Fermi liquid for small values 
of the interaction U and an insulator for large U [po[ . 
Georges and Krauth identified this transition by observ- 
ing the behavior of the double occupancy which is given 
by (D) — (n|n|), where n„ is the density of spin a, and 
of the single-particle Green's function G(r). (D)(T) dis- 
plays a minimum in the Fermi liquid regime. But when 
the system is an insulator for U > U c , (D) always in- 
creases with T. We observe the same qualitative behavior 
as in H when iV c = 1. 

We will perform the same kind of analysis to study 
the effect of short-range AFM correlations that are ab- 



sent in the DMFA. The behavior of (D) as function of 
N c depends strongly upon U. In the strong coupling 
regime, (D)(T) is essentially independent of N c (c.f in- 
set to Fig. Dh. This was expected, since for large U, 
local fluctuations dominate; they are already captured 
in N c = 1. In the weak coupling regime, the minimum 
found for N c — 1 flattens progressively as N c increases 
from 8 to 16. At N c = 36, there is a downturn in (D) at 
low T. By localizing the moments at low temperatures, 
the system gains free energy by taking advantage of the 
short range magnetic order. For the value U = 1, we are 
already deep in the weak-coupling regime of the model. 
Hence, we expect this behavior to be generic in this re- 
gion. This finding is similar to a recent finite system 
QMC simulation |£l). 
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FIG. 2. The DOS p(u>) at /3 = 32 and U = l. 

The density of states (DOS) p(u>) (shown in Fig. ||) con- 
firms the destruction of the Fermi liquid quasi-particle 
peak by short-range AFM correlations. With increasing 
N c , the gap opens fully, and the Hubbard side bands 
become more pronounced, consistent with the suppres- 
sion of (D) shown in Fig. |]. Since we work within the 
paramagnetic equations, the opening of the gap in the 
weak-coupling regime cannot be attributed to the Slater 
mechanism. 

In Fig. H we show a tentative phase diagram obtained 
with N c = 1,8,16,64. The boundary line between the 
metallic and the insulating phases is defined as the point 
where a true gap opens to the system within our numeri- 
cal accuracy; the criterion for the gap opening is roughly 
p(to = 0) < 1 x 10~ 2 . Hence for temperatures just above 
this line, p(u) displays a pseudogap. As expected, for 
U > W, the boundary line is almost independent of the 
cluster size. For U < W however, the N c — 1 line ends 
at U = U c an W. For U < U c , a Fermi liquid is stable at 
low T, For larger clusters, we don't find a quasiparticle 
peak, except when U < U c , and U c seems to approach 
zero as N c increases. This seems to indicate that there is 
no weak- coupling fixed point in the two dimensional Hub- 
bard model at half-filling. However, since an exponential 
behavior of the gap is expected in this limit, the numer- 
ical cost becomes too high at very small couplings; the 
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smallest value of U we were able to observe a pseudogap 
is U = 0.5. For instance, it was necessary that N c > 36 
in order to see a gap at U — 1 (as shown in Fig. |J). We 
believe that this gap is due to a synergism between short 
ranged order and moment formation as evidenced by the 
suppression of (D) (Fig. |l) and the enhancement of the 
Hubbard side bands (Fig. |) seen when the gap opens. 
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FIG. 3. The temperature T g the gap opens in p(u>) 
N c = 1, 8, 16, 64 as function of U. In the inset the mean-field 
transition temperature T* and the gap A as function of U, 
the linear (U > W) and exponential (U < W) fits for A serve 
as guide to eyes. 

It has been argued that the SDW picture might capture 
the description of the insulating state of the cuprates in 
the weak to intermediate coupling regime |^|| . We show 
in the inset to Fig. || the gap A and the mean-field tran- 
sition temperature, which to a very good approximation 
is also the pseudogap temperature T* characterized by a 
peak in the bulk susceptibility and the onset of the pseu- 
dogap in the DOS. T* defines the temperature where 
short-ranged AF correlations emerge. If the gap were 
due only to these correlations, then we expect A w T*. 
For example, in the SDW approximation A/T* = 1.76 
p3| . We find that A/T* strongly deviates from this pre- 
diction. In fact, we find that A > T* > T g hence, SDW 
ordering is not the mechanism that leads to the gap for- 
mation. 

We now show that the conclusions drawn from finite 
clusters are valid in the thermodynamic limit. For this, 
we exploit the complementarity between the DCA and 
FSS jl4],[l3|]. In Fig. ^, we show the imaginary-time 
Green's function G(r) at the Fermi point X = (it, 0). 
This quantity has a more rapid decay from it maxi- 
mum at G(/3/2) when the effects of the correlations are 
stronger so that the gap is more pronounced. In finite 
systems, the decay is sharper for smaller lattices while in 
the DCA the opposite occurs. This behavior marks the 
fundamental difference between the FSS and the DCA. 
At low temperatures, in FSS, the correlation length is 
greater than the lattice size. Thus the effect of correla- 
tions are overestimated for smaller clusters because these 



systems are artificially closer to criticality than a system 
in the thermodynamic limit. This tendency is reduced 
by increasing the cluster size. The situation is radically 
different in the DCA where the system is already in the 
thermodynamic limit. The DCA restricts correlations to 
within the cluster length. As the cluster size increases, 
long range correlations are progressively included. Thus, 
the effects of the correlations increase with cluster size. 
Since the two techniques must become identical in the 
thermodynamic limit, the G(t) curve in the this limit is 
bracketed by the FSS and DCA curves. 
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FIG. 4. The imaginary-time Green's function at the Fermi 
point X = (n,0) calculated with finite size QMC (filled sym- 
bols) and with the DCA (empty symbols) when U = 1.1, 
(3 = 16. 

This complementarity is also seen in the spectra. In the 
weak-coupling regime, the FSS and the DCA start from 
two different physical limits. A finite system is always 
an insulator as soon as T is less than the energy separa- 
tion between the ground state and the first excited state. 
Hence finite systems always show gaps or pseudogaps if 
T is sufficiently low. On the contrary, for N c = 1, the 
DCA is identical to the DMFA which at weak couplings 
yields a local Fermi metal. As the cluster size increases, 
one expects the two techniques, that become identical in 
the thermodynamic limit, to converge towards this limit 
from complementary directions. In the strong coupling 
regime, the DMFA has a gap too. But this it is always 
smaller than the one found in the FSS for the same set 
of parameters. Hence the convergence towards the ther- 
modynamic limit in the DCA will be from smaller gaps 
to larger, while the converse occurs in FSS. The DOS 
shown in Fig. ^ supports these conclusions. The finite 
size gap in the FSS decreases when the cluster size goes 
from N c = 16 to 64. In the DCA, there is a pseudogap 
for N c — 16 that turns into a true gap when the cluster 
size is increased to 64. Since by construction the DCA 
underestimates the gap, we can affirm that at this tem- 
perature, the gap exists in the thermodynamic limit. Its 
actual value is bracketed by the FSS and the DCA. 
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FIG. 5. The DOS p(u>) from finite size QMC (top) and 
from DCA (bottom) at U = 1, (3 = 32 and A c = 16, 64. 

Conclusion We have presented well-controlled results 
on the MIT in the 2D Hubbard model obtained by the 
DCA. We find that the double occupancy and the single- 
particle DOS show evidence of a gap for both small and 
large values of U at finite T. Since, a finite-temperature 
gap persists well into the weak coupling regime, and there 
is no long range order until T = 0, the Slater mechanism 
is likely not responsible for the metal-insulator transi- 
tion in the two dimensional Hubbard model. Instead, 
the gap is due to a synergism between short ranged order 
and moment formation, as evidenced by the concomitant 
opening of the gap and the suppression of the double 
occupancy. Since the charge fluctuations are suppressed 
at finite T by the MIT and the spins do not order until 
T=0, spin-charge separation also persists into the weak 
coupling regime. As the DCA systematically underesti- 
mates the gap formation, these conclusions are valid in 
the thermodynamic limit. The resulting phase diagram is 
consistent with Anderson's view that the effective Hamil- 
tonian for the 2D Hubbard model at half-filling for all 
U > and A ^> T is the 2D Heisenberg Hamiltonian 
|^| . Finally, our results put some constraints on theo- 
ries of high temperature superconductivity. Since the- 
ories based on weak-coupling expansions are unable to 
describe these results, their prediction in the doped case 
are questionable. 
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